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Abstract 
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Overview 


» History of investigation of mutual coupling 
» From metal horns to metahorns 


» Finite array analysis 
s: coupling between elementary sources 
o basic factors in coupling 
> element-by-element approach 
s: coupling of metahorns 
» Examples 
» Conclusion 
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Arrays of Metal Horns 
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The box horn. (a) Original rendering. (b) Array of box horns. 
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Some Other Metal Horns Used in Arrays 





Combination of box horn and Diel-core horns 
dielectric loaded "slab" horn for dual 
band operation 
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Antenna Arrays of Metal Horns 
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Metamaterials 


» Consist of structures that have combined 
electromagnetic properties as if they consisted of 
a homogeneous material with defined constitutive 
properties such permittivity (e), permeability (u) 
and conductivity (o). 
Often the constitutive properties are significantly 
different from natural materials, such as 
o permittivity (e) can be <1 
o permeability (u) can be <0 
» They are often nano- periodic structures or fraction 
of a wavelength 
> A simple and lightweight example is a wire mesh 
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X Aube 4. 
Mutual Coupolina in Arravs of Metahorns  unw 


b d 





Metamaterials 


» Can be created many different ways with a variety of 
bandwidth properties. Assume here both positive or 
negative permittivity or permeability 





Metamaterial , l4 — Physical 


u>0, <0 u>0, >0 
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Why Metahorns? 


^ Meta-horns are here defined as a hybrid mode horns with soft or hard wall 
liners. They can also be conventional metal horns containing meta-lenses. 


> Could use metamaterials which are built into the design of meta-horns. 
^ Potentially lighter weight than metal horns. 


^ Low cross-polarization possible with both soft and hard metahorns 


° Soft surface - amplitude of the E-field in both polarizations is zero at the surface eg. 
conventional transverse corrugations 


^ Hard surface - maximum value of the E-field at the surface for any polarization eg. 
longitudinal corrugations 


> Wide bandwidth is possible. 
^ Choice of metasurface permittivity versus frequency can lead to Drude-like dispersion on the 
trailing end of the Drude resonance, thereby facilitating large bandwidth. 
° This was first demonstrated on a single polarization square horn with metasurfaces of bent wires 
placed on the E-plane walls, showing over 100% bandwidth. 


^ A dual polarized square meta-horn was created with a printed circuit board metasurface that gave 
showed close to an octave pattern bandwidth. 
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Meta-horns 





C-band horn with low index 
metaliner made out of metal 
wires without dielectric on the 
E-plane walls for linear 
polarization. 


Soft Ku-band horn with low 
index metaliner implemented as 
printed circuit boards 
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Horns with Meta-lenses 


Metahorn 





Transformation optics 
Metamaterial 


Transformation optics to alter the lens 


aperture field of a horn. 
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Metamaterial Lenses in Horns 


Sidelobe Reduction Linear-to- Circular Polarizer 


-60 
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Ref: X. Chen et al., EL, 2016 


Frequency(GHz) 


Ref: C. Lin et al., AWPL, 2018 
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Mutual Coupling in Arrays 


> 


Mutual coupling is often neglected and can be useful to give 
an initial impression of the radiation pattern 

Coupling effects: 

© Sidelobes 

© Gain 

© Input match 

© Cross-polarization 

Aim here is to describe some more effects of mutual 
coupling in arrays 

By iE some of these effects, show that mutual 
coupling should be included in the original design and 
synthesis of these arrays, which corresponds to the 
selection of the complex amplitudes {A} 

Sometimes, better performance can be obtained with 
coupling included 

Applications vary from array design for radar & sensing, 
communication systems eg 5G or radio astronomy 
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Approaches 


» Floguet X 
| Element-by-element 





» Minimum scattering element X 
» Periodic structure X 
» Full array 
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Relationship of Mutual Coupling Methods 











Antenna arrays 





Infinite arrays 
^ a . 
E Periodic methods 
b 


x 
pre a Integral equation methods 


y Asymptotic methods 
Asymptotic Periodicray 
field solutions structure methods 


Y Mode matching 
Å 
4 


| 
I Á 
j 
| 








Minimum 
scattering methods 
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Scattering matrix Impedance/admittance Fourier series 
formulation matrix formulation formulation 
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PTD jet Ves 
Sa eo. | T Y k E" 
Avril 1922 —Ó" (ii IUE. Gasanbieter DN 
ST. 


RAB 
Sur l'origine de la résistance de rayonnement 


Par M. Léon BRILLOUIN 
Ingenieur eil å la Société in féependante de Téleg iphte ta ^ 
I Exposé E Tmentaire lique. Une telle pr nation n'est valable que 
Hur des courants à va r 


> 







Brillouin (1922) detailed a method of M AT 
analysis used for calculating feces || 
radiation resistance. ee MEA NN 
The P oynti ng ve ctor method was ENT ce 
described by Bontsch-Bruewitsch 
(1926) - now the textbook method 
for calculating radiation resistance - 
was used to obtain the radiation 
resistance of arrays. 





A systematized approach, called the emf method, was 
applied by Pistolkors (1929) to find the radiation resistance 
of various dipole array configurations. 


Both emf and radiation resistance methods were shown to 


be identical by means of Gauss' law by Bechmann (1931). 


Refs: Brillouin, Radioelectricite, Tome Ill, Apr.1922, pp. 147-153 
Bontsch-Bruewitsch, Annal. Physik, B.81, No. 21,1926, pp.425-453 
istolkors, Proc. IRE, Vol. 17, No. 3, 1929, pp. 562-579 
meme. Proc. IRE, Vol. 19, 1931, pp. 1471-1480 
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History - 2 


ww 


ww 


ww 
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The emf method gained wide acceptance after the 
classic paper of Carter (1932) who used it and 
reciprocity to obtain complex expressions for self 
and mutual impedances for a variety of two-dipole 
arrangements. The emf method became the basis of 
the element-by-element approach 


Floquet method was adopted from solid state 
physics in the 1950s and 60s. 


Periodic structure method Wasylkiwskyj and Kahn in 
1970 


Full array analysis became possible with 
development of numerical methods and computer 
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er in the 1980s & 90s. 


: ahn, IEEE TAP, AP- 18, Nov. 1970, pp. 741-752 


Mutual Coupling in Arrays of Metahorns 





Element-by-Element 


P 


Ww 


wv 





In this approach, an impedance or admittance 
matrix is created by considering in turn each 
element of the array in the presence of another 
element, with all other elements of the array 
removed either by physical means or by open or 
short circuiting the terminals. 

For wire antennas all except the two being 
considered are left open circuited and the mutual 
impedance is calculated. 

For slot and waveguide antennas the mutual 
admittance of the two elements is obtained with 
the remaining elements short circuited. 
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Simplified Models of Electric & 
Magnetic Dipoles 









[rtt Me) 
[E oon 


(a) 


jk e ani 
dE, = -17 —u— dm*R 


dg = À R x dEn 





n=wave impedance of the medium -jë 
ay 


Ho 
(b) To - jé =wave impedance of free space 
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Elementary Sources 


dp — J.4 5" 
dm = M.dS" 


Electric current element 
jk o JKR E ^ 
dH, - 7| —,— dp KR dE, —7, dH, xA 


Magnetic current element 


jk e ^5 b 1 b 
dE = — —_ d a x R dM = — dE R 
m 4r R a 7» m ^ 


Total field from both sources 





dE - dE, --dE,, dH - dH, -dH,, 
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Extension to Two Elemental 
Dipoles - Mutual Impedance 






Z, Z4 Z5 
g. 
uS X) 
l2 
£paT-7-7--------------7 QA 
Ms Dipole 2 
PN 
| 
I1 | 
| | 
7 | 
f X, X4 
^ Xo2 
Dipole 1 


From the theory for short dipole, the fields due to 1 are: 






Lidli onla i . ] 
dE, = lo ^1*^1 oim 0, sin Ó, Eat tf sont. dH, — — f, xdE, 
I 4r n n* n* 7] 
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O 


Elemental Dipole Coupling -2 


The emf across dipole 2 due to 1 is 


The mutual impedance is defined 
dV,, dE,, :z, d^, 
£4, = ooo 
L L 





Similarly for dipole 1 due to 2 


dV dE ,, :Z, df, 
fp o 
2 2 


From Lorentz reciprocity theorem 
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and therefore Zp = Z2 





Elemental Dipole Coupling - 4 


Mutual impedance of the dipoles: 


Medes s, 4 jk 
Zz = 
47 


l 
sin 6, sin(@,, —a@)+— (cos o: (1—3cos^ 0, ) - sin a sin 0, cos 6, (2 cos 6., — ») 
a 52] 


E 

Parallel dipoles & H-plane (a=0 ): | 
Z, = n, dt d£ » okt jk D 
47 1 Vi 


In echelon arrangement (ie. end-to-end and in E- plane): 


— Fo ae? e E 


Z7 
2m 27 


Booker's relation: 


(Superscript efrefers to electric field) 


25 
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Basic Facts About Coupling 


» The essence of mutual coupling 
o° source 
s: dimensions of transmitter & receiver 
o distance 
> polarization 
o orientation of receiver 
» From these, by DET all mutual coupling 
can be interpreted [V] = [ZIV] or Ul - ITIIV] 
» In terms of voltage and current waves 
IV] = [al + [b] and [= [a]- [b] where [2] = [s][a] and 
[S] = [1 +Z] iz -1] [S]= 5 +Y] >- Y] 
» The coupling coefficients are S; where i+} 
> E example, if Yi—^1 /n, and Y; « « Yi then S;z:- n, 
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Element-by-Element Approach - Dipoles 


A A V NE 


Antenna #1 Antenna #2 Antenna #n 


J; =current on #1 in the absence of #2 

E,, =electric field in the absence of #2 due to #1 
J, = induced current at #2 

Any other elements are open circuited 


Based on the earlier discussion of elemental dipoles, namely 


dV,, dE,, -z, dé, 
Za = rr 





I A 


summation of all such dipoles over the surface of antenna #2 results in 
autual impedance 
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Element-by-Element Approach - Slots 
vA D M; M =0 


n 





Antenna #1 Antenna #2 Antenna Zn 


M, —magnetic current on #1 in the absence of #2 
H,,=magnetic field in the absence of #2 due to #1 
M, = induced current at #2 


All other slots are short circuited 


From duality with previous case, the mutual admittance is 


Note: The expression for mutual impedance and its dual for mutual 
admittance are variational or stationary with respect to second-order 
changes in the current/voltage 
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Methods 





Array size 


» Mode matching 
» Method moments 


LL 


> Characteristic modes | [o 
» Field correlation v. v 
» Numerical methods Dimensions 





> FEM, FDTD etc 
» Hybrid method 
> MoM, mode matching and field correlation 
> MoM & characteristic modes 
> Numerical methods & mode matching 
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Mode Matching - Apertures in Ground 





T T 
I | | 
— 
a, we x 
B s^ LlFl^g | 
my” bc, 
| 
| [4 —— — s———— 
M (1) f — ; (i) i Jy e 3 tjy Z +1/2 
j n B t Jy yi z L) 2. = dr pte pi 
E = S ae Tpi? tbe a )e,;05)Y, H, = X (ae b € )h, Go, y)Y, 
p=l p=l 
/pi 7 Ppi - J*pi is the propagation constant of mode p in section i. 


jl. d$ ey xhg-2-26,, 


Ypi is the wave admittance of mode p in section | and “is the Kronecker delta. 
The external magnetic field for magnetic currents on the aperture is 


Hext =| M -G“Y(R|R’) ds 
ZUTS -~ 


Sext 
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Mode Matching - 2 


H~ 
we 
M, 


\ 


S, 


In the aperture, starting from the continuity of the magnetic field on S, 
H;* -H;" (2, - Ht" (a); j 5L... N, 
where a; are unknown coefficients in which M, is expressed 

M, =YV a, 

The continuity condition is now 
R = Sara )-H®™(a@,)|+ H™ 

Now weight this residual by a function p, and integrate over S, setting the 
result to zero ie. (5. R)- [[as' AR -o 


S 
= ffas s EEn e-ur eph ul =0;i=1,...,M, 


Special cases: 1. when Z, =6(«—-x,)6(y—y,) the approach is point-matching or 
collocation, 2. When N.=M, & {,}={a,} the approach is Galerkin's method. 
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Mode Matching in the External Region 
Outside a conductor on the (u1,u2) plane, the external region is 


represented by ` 
v» -p;]- i [4s & mia) |- ffas g È Ha) 


lec perm 
S, 5, 


E Dy. (m, n) A] Y, Y; 


where the normalized mutual admittance is expressed 








Once V has been found from a solution of 
(ye + y~ Jv 2 1" 


the mode amplitudes can be found from 


y"? N, 
(an +bn)= == 2 (a hi, and the scattering matrix 
j=l 
b=Sa where — S" -(I«y) a-y) 
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Method of Moments: Example Non- 
contacting Horns 


P P 
— at aw 
| J l| \\ | 
F-- / 7 Sit Rr piri Fa : waveguide 1 : M, 
| My 91 | 2 
\ \\ f ! n 
Weca A i U 
Z-0 
pl 
rr ERR 
| Jo I] | | | 
n M5 | 2) 22 waveguide 2 tM, 
\ \\ j) i 
CE S O E v |o 
Zo-0 P 









Ref: Cooray & Bird, Radio Science, Vol. 35, 2000, pp. 5 
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Formulation - 2 


Satisfaction of the boundary conditions gives: 


<W {nj x (By + Ej;g)»-0 on SjUA, i-1,2 


ti 


zw 


mau? 


[ñi x (Ha + Hi)]»-«W fa; x Hi» on A, i-1,2 
where the inner product is defined as: 
(A, B) = | (A - B) ds. 
J 5 
Application of weighted residuals results in: 


[z d [DU] 


[[Z = 





m 
uni je") TI iu DA 


Solve for unknowns [6] through solution of the linear equations 
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General Remarks on 
Coupling 
The Greens function can be approximated for large separations or large 


bodies. This enables approximate asymptotic formulae to be derived that 
are applicable at large distances. 





. Thus: 


Y, o (radiated field from aperture 1) m: (radiated field from aperture 2) ocu X 


[2 
(mutual admittance of magnetic dipoles in situ) 


Ol 


The normalized admittance is y, =Y," Y Y,” 
mode in aperture i. 


If the apertures are in a ground plane, then use the mutual admittance of 
magnetic dipoles in a ground plane. 


If the apertures are in a metallic cylinder, use the mutual admittance of 
magnetic dipoles on a cylinder. 


abplane coupling of two circular waveguides 
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is the wave admittance of the 






Loading Horn Walls 
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Design Techniques for Metahorns 


» The boundary conditions for the surface impedance (Z'^and Z!) of an ideal horn 
liner that are required to support the soft, hard and balanced hybrid are as follows: 





ype OT 5urrace Soundary Conditions 


Soft ane 
Hard D rmm so 
Balanced Hybrid Mode XTEXTM — — 06 
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Synthesis of Liners for Metahorns 


» Engineering the metamaterial liner of a horn antenna to possess the dispersive properties 
the suitable mode over a broad bandwidth is extremely challenging 

^ powerful synthesis procedures are used, which are typically based on global stochastic 
optimizers such as the genetic algorithm 

>» the cost functions utilized in the genetic algorithm optimization procedure to achieve 
broadband soft, hard, and balanced conditions are given below: 








cost n2 Cum em 
SO 
frequencies AM 
cost m (rs au 
u frequencies xe 








= TE yTM 2 
COStpalanced — X l (X x + no) 
frequencies 
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Example of Metamaterial Liner for a 
Horn 


Air 


Patterned Screen 


» (a) Side view of the 
metamaterial liner unit cell 
showing patterned metallic 
screen on top of a thin 
dielectric substrate with a 
metal backing and a via; 


» (b) an exploded view of a 
Single pixilated 
metamaterial unit cell 
geometry; and 


» (c) isometric view of 
optimized metamaterial. 





Ref: Q. Wu, et al, IEEE TAP, vol. 60, pp. 3594 -3604, 2012 
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Rectangular Meta-horns 





Trifurcated horn Slab-loaded H- plane Slab-loaded E-plane 


Ref: Chan et al. 1980 Ref: Tsandoulas & Fitzgerald 1972 Ref: Lier & Shaw 2008 
- Improved efficiency - Improved beam symmetry 
- Reduced BW in horizontal plane — Reduced sidelobes in E-plan 
- Higher overall mutual coupling - Lower overall mutual couplit 
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C-band Horn with Wire Grid Liner 


120 
| 100 
AN 
| 60 
40) 
20 


Aperture distribution 5.85GHz Co- pol (solid), cross-pol. (dashed) at 5.85GHz 


E field magnitude Vini 
Gan (dB) 





à (degrees) 
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Slab Dielectric Horn Array 


um 
P 


2.5 





+ cs 
B 
s ze ` 
| g 
E FE 
c 3 
V H 
X ze 
& |. 
RS : 
ab s 
F a A. 2 1 g 1 2 
Asimuth (degrees) 
3 Directivity contours for New Zealand transmit 
beam. Contour levels are in 5dB increments 


up to 35dBi, then 1dB thereafter. 





Radiation patterns of slab loaded 
horn in 45 deg plane at 12.5GHz. 


reflection coefficient, dB 





b A 10 Fa 14 


frequency , GHz — 


Reflection coefficient as a function of frequency 

for dielectric loaded waveguide a=b=22.86mm, 

d=3.46, ¢,=3.84, ¢,=1. 
---— - theory (LSE,, only} 


Ro a Serie nl 





theory (15 mates) 
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Principal Plane Patterns 


E-plane H-plane 


Power (dH) 
“15 e 





0 30 60 90 
Angle (deg) 


Theory solid curve; Expt. dashed 
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E-plane Loaded Horn - Aperture 
Distribution 


Time harmonic solution with implied factor exp( jat— j7z) 
A TM, (hybrid) mode with n>0: 


, 2 
— nz [n 
EV = L C, cos(k,x)sin 1 
UA M OEE, b b 
N , 2 
— nz a .[n 
E EL up m C,cos k, —-—x ||sin 2L 
WEE ,-9 b 2 b 
nz n 
———— kci in s}eos{ 22| 
b 


es Pee 


k£» sin kw 
k; 2 Jo eus - Y. and C, --C m ; 


1 
k€ sin kzt 








Note, if ¢.,->e¢,, E, & E, ->0 
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Radiation from Loaded Horn 


y is the propagation constant. The dispersion relation between y and frequency k-ovc for 
LSM; mode is obtained from 


hz hı 
— tanh, t = ——tanh,w 
£y2 Er1 


Propagation constant is approximately given by 


ye efisen -02-i ns WH) (2) 


The far-fields are in the same form as air-filled pyramidal horn ie. 





jk e 


E,(r,0,Q) BE 
Am r 


N (0,9) (1-- cos 0) cos ø 


— jkr 


i | 
Eg 10.9) =- : : N .(0,9) (1- cosO)sin where 








a/2 
—j nm? 2nny'A | C ! a 
N,(u, v) = — Jp + (=) | dy! ei(2nvy'-ky""/2pH) cos ajl dx'ei(rux-kx:*/2pE) cos L (5 = x) 
(0£g b » a 2 
—-a — 


Ww b/2 
$ C dx'ei(2Ttux'-kx1?/2ppg) cos(A x^) 4 C2 | dx! ej (2mux'—kx!? /2pz) GS L (= B x) 
Er1 as €r2 n 2 






1. 1. " 
u=—sindcosp v=—sinOsing & Pen are radii of curvature 
\ in E- and H-planes of horn 
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Aperture Taper - Metamaterial Lining 


Consider the case of a horn with aperture a=b= 6.1 cm, p;=33.1 cm 
if ¢ ,=0.45 and up.=1 and t=1 mm, the E-plane aperture distribution is: 


E-plane pattern 


Aperture 
2.0 
Y 15 
= 
a = 
E 10 S 
D 
z 
£ 


` -1.0 -0.5 0 05 1.0 
X, (A) 





0. 20. 40. 60. 80. 100. 
Angle, (deg.) 
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Standard Gain Pyramidal Horn with 
Low Sidelobes 


—@— airhorn 
—h— metahorn 


T" 4— airhorn at 12.9GHz 
—*— airhorn at 13.2GHz i M 


relative power level, dB 
relative power level, dB 





-90 -45 0 45 90 -90 -45 0 45 90 
angle, deg angle, deg 


Frequency 12.9 GHz 
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Meta-horns in a Ground Plane 












Transformation optics region 
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Analysis of Mutual Coupling 








ds'M,(R’). | | dSGV? (R|R;) M;(R) 
S; Si 


where 


G (RIR) = a x res d ve) G,(|R — R'| 





RR |= Vex? + 0-¥)P + (@-2 






M(R) = 27 X Eqy(R) 
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Coupling in Rectangular Horns 


Aperture fields with metamaterial loading 


Cosine-shaped Aperture Fields 
either E- or pen 
Uniform Aperture Fields -2 
either E- or H-plane 


aie O05 A, 


EB 


SA 





- AU 





- TE 40777 TE 45 

in c ` 
T H S 
t & | ` 
a a 
5 E 1 Ar? 
= E 
a C 
i 

m En 
c i 
2 4 
& a 
o 
c -5 r 

0-6 Q7 D- A 0-9 1-0 1:1 1:2 123 
"4[X) 
Fok -5j ——————a--e-4--.- I i i 
0-65 Oh? -Å 0-8 t: 1*1 Ie? LX 


Asymptotic decay is 1 /r Asymptotic decay is 1/r? 
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Conclusions - ] 


» Mutual coupling is a necessary consequence of 
any antennas in close proximity 

» Mutual impedance/admittance is a summation 
of distributed sources from one antenna on 

another with appropriate Green's fun ions ie. 






» Approximate models for mutual coupling in 
horns are useful in the initial design process. 
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Conclusions - 2 


» Where possible modern computational and/or 


ww 


{v 





measurement tools should be used to include 
mutual coupling in detailed design 

Design for mutual coupling from the outset will 
save time in the long run, maximize performance 
and remove uncertainties traditionally associated 
with mutual coupling effects eg. array blindness. 
Mutual coupling between metahorns can be 
handled using similar techniques as metal horns 
but using appropriate aperture currents. 


o Derived by field transformation or modes of loaded 
structure 


> Asymptotic decay depends on taper of aperture fields. 
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